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ABSTRACT 

NGC 6369 is a double-shell planetary nebula (PN) consisting of a bright annu- 
lar inner shell with faint bipolar extensions and a filamentary envelope. We have used 
ground- and space-based narrow-band optical and near-IR images, broad-band mid-IR 
images, optical long-slit echelle spectra, and mid-IR spectra to investigate its physical 
structure. These observations indicate that the inner shell of NGC 6369 can be de- 
scribed as a barrel-like structure shape with polar bubble-like protrusions, and reveal 
evidence for H2 and strong polycyclic aromatic hydrocarbons (PAHs) emission from a 
photo-dissociative region (PDR) with molecular inclusions located outside the bright 
inner shell. High-resolution HST narrow-band images reveal an intricate excitation 
structure of the inner shell and a system of "comet ary" knots. The knotty appearance 
of the envelope, the lack of kinematical evidence for shell expansion and the apparent 
presence of emission from ionized material outside the PDR makes us suggest that the 
envelope of NGC 6369 is not a real shell, but a flattened structure at its equatorial 
regions. We report the discovery of irregular knots and blobs of diffuse emission in low- 
excitation and molecular line emission that are located up to 80" from the central star, 
well outside the main nebular shells. We also show that the filaments associated to the 
polar protrusions have spatial extents consistent with post-shock cooling regimes, and 
likely represent regions of interaction of these structures with surrounding material. 
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1 INTRODUCTION 

The planetary nebula (PN) NGC 6369 possesses a double- 
shell morphology consisting of a round bright inner 
shell and a fainter filamentary o uter shell or envelope 
(|Schwarz, Corradi, &; Melnickl [l992). Its most remarkable 
morphological features are two disparate extensions of the 
inner shell along the east and west directions; the latter 
appears as a large, filamentary blister or ansae, whilst the 
former is a bifurcated structure. 

Using low-dispersion spectroscopy, iMonteiro et al.1 
<|2004h modelled the peculiar morphology of NGC 6369 
as a tilted collar of gas and two symmetrically located 
caps of emission. This interpretation is consistent with the 
spheroidal shell and faint bipolar extensions proposed by 
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ISteffen fe Lopezl (2006) based on a single long-slit high- 
dispersion spectrum obtained along the nebular major axis. 
These simplified models do not deal with the bifurcated 
structure to the eastern side of the inner shell, nor they 
consider the filamentary appearance of the envelope. These 
have been propos ed to represent fast low ionisation emission 
regions fFLIERs. iHaiian et al. 1997), although evidence for 
fast expansion velocities is lacking. 

The H 1 and He 11 Zanstra temperatures of NGC 6369 
have been reported to be T Z (H 1) ^ T z (He 11) ^ 70,000 K 
(e.g., lPhillipsl i2003), comparable to estimates deriving from 
the energy balance method ([Pottasch &; Bernard- Salas 
2008). The similarity of the H I and He II Zanstra tem- 
peratures suggests that the nebula is optically thick to H 
ionizing radiation, and thus that significant amounts of 
molecular material and dust may be present outside the 
ionisation bound optical nebula. This is consistent with 
the ISO (Infrared Space Observatory) indications of emis- 
sion from polycyclic aromatic hydrocarbons (PAHs) bands 
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Figure 1. {left panels) EST WFPC2-WF3 images of the inner shell of NGC6369 in the Ha, [O III] A5007, and [N n] A6583 emission 
lines, and (right panels) [N n]/Ha, [N n]/[0 III], and [O m]/Ha ratio maps. 
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(jCohen fc Bar low 2005) that would be produced at a photo- 
dissociation region (PDR) located outside of the ionized 
shell. There is also evidence of dust emission from NGC 6369. 
ISO spectroscopic observations of NGC 6369 have showed 
it to be one of the few evolved PNe to possess the 21 a nd 30 
/im emission features (iHonv. Waters. &; Tielensll200lh that 
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spectral features is unclear, although it has been suggested 
that the 21 fim emission may arise from an extraordinarily 
wide range of agents, including hydrogenated fullerenes, 
nano diamonds, SiC+Si02 g rains, and TiC nanoclusters (see 
e.g.. lZhang. Jiang. &; Lill2009l . and references therein). The 
carbon dust origin in NGC 6369 is consiste n t with its high 
C/O ratio, -2 (jZuckerman fe Allerl Il98& lAller fe Kevesl 
1 19871 ). 

In order to investigate the kinematics and three- 
dimensional structure of the main nebula and envelope of 
NGC 6369 and to confirm the occurrence of an outer enve- 
lope of molecular material, we have acquired new narrow- 
band near-IR H2 and optical images, and long slit high- 
dispersion spectroscopic observations that have been exam- 
ined in conjunction with archival mid-IR ISO spectra and 
Spitzer images and spectra, as well as optical narrow-band 
HST images. The observations and archival data are pre- 
sented in 321 and the results on the morphology, and mid-IR 
spectroscopy are described in 33 and 311 respectively. The 
kinematics and the discussion are presented in 33 and final 
conclusions are presented in Jj6j 



2 OBSERVATIONS 
2.1 Optical imaging 

Narrow-band images in the [N 11], Ha, and [O ill] emission 
lines were obtained on September 1, 2008 using ALFOSC 
(Andalucia Faint Object Spectrograph and Camera) at the 
2.56-m Nordic Optical Telescope (NOT) in the Observatorio 
de Roque de los Muchachos (ORM, La Palma, Spain). The 
detector was a 2048x2048 EEV CCD with plate scale 0719 
arcsec pix -1 , resulting in a field of view (FoV) of 6'5x6.'5. 
The central wavelengths and bandwidths of the filters are 
A C =6589A and AA=9 A for [N 11], A C =6568A and AA=8 A 
for Ha, and A C =5007A and AA=8A for [O III]. The spatial 
resolution achieved during the observations, as determined 
from the FWHM of stars in the FoV, was tf/95. Either two 
or three individual frames with integration times of 300 s 
were taken for each filter, leading to total exposure times of 
900 s for the [N 11] image, and 600 s for the Ha and [O ill] 
images. The data were bias-subtracted and flat-fielded with 
twilight flats using standard IRAF0 routines. 

Narrow-band HST WFPC2 images of NGC 6369 in the 
[N 11] A6583, Ha, and [O ill] A5007 emission lines were 
retrieved from MAST, the Multimission Archive at the 



1 IRAF, the Image Reduction and Analysis Facility, is dis- 
tributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy (AURA) under cooperative agreement with the National 
Science Foundation. 



Space Telescope Science Institut43 (Prop. ID: 9582, Hub- 
ble Heritage Observations of NGC 6369, PI: Keith Noll). 
The nebula was registered by the WF3 detector with a 
pixel scale of 0V09. Four individual exposures were ob- 
tained using the F658N (pivot wavelength A P =6590.8A, 
AA=28.5 A), F656N (A p =6563.8 A, AA=21.5 A), and F502N 
(A P =5012.4A, AA=26.9A) filters, resulting in total expo- 
sure times of 1,040 s for the [N 11] A6583 image, and 640 s for 
the Ha and [O ill] A5007 images (Figured). 



2.2 Near-IR imaging 

Narrow-band H2, Br7 and continuum K c near-IR images of 
NGC 6369 were obtained on 2010 June 27 using LIRIS, the 
Long- Slit Intermediate R esolution Infrared Spectrograph 
([Acosta-Pulido et al.ll2003h , at the Cassegrain focus of the 
4.2-m William Herschel Telescope ( WHT) on Roque de Los 
Muchachos Observatory (ORM, La Palma, Spain). The de- 
tector was a lkxlk HAWAII array sensitive in the spectral 
range from 0.8 to 2.5 fim. The plate scale is 0'/25 pix -1 and 
the FoV 4. / 27x4. / 27. The narrow-band filters have A c =2.1218 
fim and AA=0.032 fim for H 2 (1-0) S(l), A c =2.1658 fim and 
AA=0.032 fim for Br7, and A c =2.270 fim and AA=0.034 fim 
for K c . 

Ten 60s exposures were taken for each filter, for total 
effective exposure times of 600 s. The jitter option was used 
to increase the sampling of the sky: the telescope point- 
ing of each exposure is shifted by a few pixels, thus raster- 
ing the nebula to different locations on the detector. Each 
series of observations on the object was followed by simi- 
lar series of observations on adjacent blank sky positions. 
The reduction of the LIRIS data was carried out using 
the dedicated software LIRISDR (LIRIS Data Reduction 
package), a pipeline for the automatic reduction of near-IR 
data developed within the IRAF environment. The reduc- 
tion performed by LIRISDR includes standard and addi- 
tional non-standard steps such as bad pixel mapping, cross- 
talk correction, flat-fielding, sky subtraction, removal of re- 
set anomaly effect, field distortion correction, and final im- 
age co-addition. The continuum contribution to the H2 and 
Br7 emission was subsequently removed using the K c im- 
age. In this process, we did not find necessary to scale the 
images since the H2, Br7, and K c filters have very simi- 
lar bandwidths and transmissions. The spatial resolution on 
the WHT images, as determined from the FWHM of stars 
in the field of view, is ^O'/S. The narrow-band continuum- 
subtracted H2 and Br7 images of NGC 6369 are presented 
in Figure [2] 



2.3 Mid-IR imaging 

Spitzer Infrared Array Camera (IRAC; iFazio et all [2004) 
images of NGC 6369 (Figure 3) were retrieved from the 
NASA/IPAC Infrared Science Archive (IRS A). The images, 
obtained on 2005 September 22 as part of Program 20119 
(The Darkest Cloud, An IRAC/MIPS Survey of the Pipe 
Nebula, PI: C.J. Lada), have been processed as described 
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Figure 2. WHT LIRIS H2 (Ze/it) and Br7 (right) background-subtracted images of NGC6369. The dark and light diffraction spikes are 
a result of the rotation of the focal reducer during the acquisition of the images. The envelope and inner shell are labeled. The arrows 
point to knots and blobs of diffuse H2 emission. In the envelope, five knots, marked as knots Kl to K5, have been highlighted. 



in the IRAC data handbook The observations employed 
niters having isophotal wavelengths (and bandwidths AA) 
of 3.550 /im (AA=0.75 /im), 4.493 /im (AA=1.9015 /mi), 
5.731 /im (AA = 1.425 /im), and 7.872 /mi (AA = 2.905 
/ im) . IRAC spatia l resolution varies between ~l"7 and ~2"0 
( Faz io et al.l [2004), and is reasonably similar in all IRAC 
bands, although the diffraction halo at 8 /mi is stronger than 
at the other bands. 

2.4 Long-slit high-dispersion optical spectroscopy 

Long-slit high dispersion optical spectroscopy of NGC 6369 
was obtained on July 17-1 9 2008 using the M anchester 
Echelle Spectrometer fMES. lMeaburn et alJl2003h mounted 
on the 2.1m (f/7.5) telescope at the Observatorio As- 
tronomico Nacional de San Pedro Martir (OAN-SPM, Mex- 
ico). A lkxlk SITe CCD was used as detector with a 2x2 
on-chip binning, resulting on a plate scale of O'.^pix -1 and 
a dispersion of 0.1 A pix -1 . Since MES has no cross disper- 
sion, a AA=90A bandwidth filter was used to isolate the 
87 th order covering the spectral range which includes the 
Ha and [N 11] A6583 lines. At this spectral order, the slit 
width of 150- /im (2") that was set during the observations 
corresponds to a spectral resolution of c_ 12kms _1 . 

Five long-slit spectra, labeled from 1 to 5 on Figure [H 
were obtained to map the kinematics at different regions of 
the nebula. The position angles (PAs) for slits 1 to 5 are 
+74°, +80°, -22°, +22°, and -86°, respectively. The expo- 
sure time was 1,200 s for the spectrum at PA=+80° (slit # 
2) and 1,800 s for the rest. The seeing during the observa- 
tions, as determined from the FWHM of stars in the FoV, 
varied from l'!2 to 2'/0. The spectra were wavelength cali- 
brated with a Th-Ar arc lamp to an accuracy of zblkms -1 . 

2.5 Mid-IR spectroscopy 

Spitzer Infrared Spectrograph (IRS; iHouck et al.l [2004) 
mid-IR spectroscopic observations of NGC 6369 obtained 
through the Spitzer Program 40115 (Dual Dust Chemistry 

3 http://ssc.spitzer.caltech.edu/irac/dh/iracdata_handbook_ 
3.0.pdf. 



in Wolf- Ray et Planetary Nebulae, PI: G. Fazio) were re- 
trieved from IRS A. The observations consisted of short 
(5.1-8.5 /im) and long (7.4-14.2 /im) wavelength spectra 
taken with the Spitzer IRS Short-Low module 2 (SL2) 
and 1 (SL1), respectively. Both modules have rectangular 
apertures ^3 / /6x57 // that, for these particular observations, 
were oriented along PA=177.18°. The observations were ac- 
quired at closely similar positions a(J2000.0) = 17 h 29 m 20*37 
and <5(J2000.0)=-23° 45' 24'.'55 (positions A & C), and 
a(J2000.0) = 17 h 29 m 20! 44 and £(J2000.0)=-23° 45' 43 / /55 
(positions B & D), resulting in overlapping slits covering 
the nebula along its minor axis. The spectral resolution of 
the observations, A/AA, varied between 64 and 128. The 
Spitzer SL1 & SL2 spectra of the source are illustrated in 
Figure [9] where we show the results for overlapping aper- 
tures across the centre of the source (labelled A to D). All 
positions cover the inner shell, but positions A and C map 
the north portion of the envelope and positions B and D the 
south portion. 

Similarly, we mad e use of mid-IR ISCE Short Wave 
Spectrometer (SWS; Ide Graauw et al.| [l996) spectroscopic 
observations of NGC 6369 acquired on 1997 February 14 as 
part of program 45601901 (PI: S. Gorny). ISO SWS provides 
spectra in the short (2.38-12.1 /im), medium (12-29 /im), 
and long (29-45.3 /im) wavelength ranges, acquired through 
apertures with sizes 14 // x20 // , 14 // x27 // , and 20 // x22 // , re- 
spectively. The ISO SWS spectra have an exposure time of 
1,140 s, and were centred on a(J2000.0) = 17 h 29 m 20!78 and 
£(J2000.0)=-23° 45' 32 / /.l. 



3 MORPHOLOGY 

An inspection of the optical and infrared images of 
NGC 6369 (Figures [1] [2] and [3} confirms the three morpho- 
logical components previously described in the main nebula: 
a bright inner shell, two extensions, and an envelope. Faint 

4 Based on observations with ISO, an ESA project with instru- 
ments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom) and 
with the participation of ISAS and NASA 
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Figure 3. Spitzer IRAC images of NGC 6369 overplotted by contours. The contours correspond to levels of the intrinsic surface 
brightness, SB, defined by SB n =Ax lOO" 1 ) c — B MJy sr" 1 , where n is the contour level (n=l corresponds to the lowest and outermost 
level), B is the background level, and A and C are scale parameters. The contour parameters [A, £>, C] are given by (1.6, 1.018, 0.199) 
at 3.6 /mi, (1.45, 0.792, 0.196) at 4.5 /mi, (7.3, 5.264, 0.17) at 5.8 /mi, and (20.5, 15.341, 0.173) at 8.0 /mi. We note that the blob of 
diffuse emission in the 8 /mi image inside the inner shell is an artifact, an echo of the brightest region of the rim. 



blobs of diffuse emission are also detected outside the main 
nebular shells. Next we describe in detail the morphology 
of the main nebular shell and this newly detected extended 
emission. 



3.1 The main nebula 

3.1.1 Optical morphology 

The HST WF3 Ra and [O III] images of the inner shell and 
envelope of NGC 6369 presented in Figure [T] show a thick 
(~5'!b) annular shell with inner radius ~8'/3. The surface 
brightness is noticeably uneven, with patches of bright and 
diminished emission. Several dark knots and lanes, reminis- 



cent of those seen in other PNe (e.g., IC 4406 and NGC 6720, 
lO'Dell et aLll2002h are also detectable. 

The [N n] image reveals a dramatically different view 
of the inner shell, with the Ha and [O in] dark knots 
and lanes turning into bright [N n] spots and filaments, 
and an intricate system of [N n] filaments filling the in- 
nermost regions of the shell. Radial filaments emanate out- 
wards from most of the [N n] knots, which are thus revealed 
as "tadpole" or "comet ary" features very similar to those 
seen, e.g., in NGC 239 2 jO'Dell et alJl2002h and NGC 7354 
([Contreras et al.ll2010h . Interestingly, there is a bright ring 
of these comet ary knots, at the location of the outer rim of 
the bright Ha and [O in] shells at a distance of ~12 / /4 of 
the central star, whose tails extend at least up to ~20". 

The ratio maps in Fig. [1] reveal that the bright rim has 
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Figure 4. Slits positions used in the long-slit echelle spec- 
troscopy of NGC6369 superimposed upon the HST WFPC2- 
WF3 [N n] image. 



high excitation, with the lowest [N n]/Ha and [N n]/[0 in] 
values at radial distances 8"- 12", but the regions within this 
shell, closer to the central star, have anomalously low ion- 
ization. At the location of the ring of knots, at ~12 / /4 from 
the central star, the "comet ary" features are highlighted in 
the [N n]/Ha and [O in] /Ha ratio maps (Fig. [Q, showing 
high ratio values in the former and low ratio values in the 
latter. The tails of these knots define an annular region with 
outer radius ~18". 

The peculiar extensions of the inner shell of NGC 6369 
([Schwarz, Corradi, &; Melnick 1992) along PA^63° are best 
seen in the Ha and [N n] images. The western extension 
is reminiscent of a lobe or a large "ansa", whilst the east- 
ern one displays a complex morphology best described as 
a bifurcated structure. The ratio maps indicate that these 
features have low excitation, with values of [N n]/[0 III] up 
to 20 times larger than those in the higher excitation regions 
of the inner shell, and up to twice those of the ring of knots 
at the rim of the inner shell. 

The envelope is barely detected in [O in], and only a 
little brighter in Ha. The [N n] emission is brighter and 
reveals it to be composed of an ensemble of knots, blobs, 
and filaments. The appearance of these knots and their tails 
is less defined than those in the inner shell. Overall, the 
emission is mostly distributed within a fragmentary annular 
region with inner and outer radii 25 // -34 // , as illustrated by 
the NOT composite-colour picture of NGC 6369 presented 
in Figure [SJ-Ze/t. 

3.1.2 Near-infrared morphology 

The morphology of the bright inner shell and envelope of 
NGC 6369 in the near-IR Br7 line is similar to that in Ha 
(Fig.^i right). The lack of differential extinction between the 
optical and near-IR lines of H I suggests little dust content in 
the nebula along the line of sight. By contrast, the image in 
the near-IR H2 line (Fig.^left) offers a very different view of 
the main shell of NGC 6369. The region interior to the bright 



inner shell, which showed anomalously low excitation emis- 
sion, is revealed to be filled with H2 emission. This region 
is then encompassed by two concentric annular regions, the 
first one where the H2 emission is very low, and the second 
one, with radius ~13", where we find bright H2 emission. 
The interior annulus is coincident with regions of high exci- 
tation in the [N 11] /Ha and [N n]/[0 ill] ratio maps, whilst 
the exterior, H2-bright annulus is coincident with the rim of 
the inner shell defined by the heads of the [N 11] "cometary" 
knots. This provides evidence that the [N 11] knots and fil- 
aments are associated with H2 emission, i.e., they include 
significant amounts of neutral material. 

The extensions of the inner shell are not detected in H2 , 
but the northwest quadrant of the envelope and some bright 
knots exhibit H2 emission. Surprisingly, the radius of the 
northwest arc detected in H2 is ^25", i.e., this H2 emission 
is interior to the annulus of [N 11] emission in the envelope. 
As for the H2 knots, some (e.g., K3 and K5) have easily 
identifiable [N 11] counterparts, while for some others (e.g., 
K2 and K4) the identification is hampered by their spatial 
coincidence with the [N 11] -bright eastern and western exten- 
sions of the main shell. Among the H2 knots detected in the 
envelope, the morphology of knots K2 and K3 is quite no- 
ticeable, with a bright head and two tails pointing outwards 
from the central star. 



3.1.3 Mid-infrared morphology 

The mid-IR morphology of NGC 6369 is illustrated by the 
"unsharp" Spitzer IRAC colour-composite picture presented 
in Figure fright. In this picture, the bright optical shell 
is the brightest feature, whereas the outer envelope and 
east and west extensions are fainter but still noticeable. De- 
spite the reduced spatial resolution, the mid-IR morphology 
of NGC 6369 is consistent with that observed in the HST 
(Fig. [Q) and NOT (FigEJ/e/fc) optical images. The redden- 
ing of the eastern and western extension of the main shell 
and outer envelope in this colour-composite IRAC picture 
implies a higher relative contribution of 8.0 fim. emission at 
these morphological components. This is further evidenced 
by the contours in Figure [3j the inner ring is intense in all 
of the IRAC images, whereas the surface brightness of the 
envelope becomes relatively brighter at longer wavelengths, 
in the 5.8 and 8.0 /im bands. 

The mid-IR surface-brightness profiles extracted from 
background-corrected IRAC images through the minor and 
major axes of the nebula (Figure [6} confirm this result. The 
shape of the surface brightness profile interior to the bright 
main shell is similar in all of the channels, with flux in- 
creasing with wavelength. By contrast, the 3.6 and 4.5 fim 
surface brightness profiles fall-off sharply at ^25", whereas 
the emission at 5.8 and 8.0 /mi show broad skirts up to 
^70". This becomes in a sharp increase of the 5.8/xm/4.5/im 
and 8.0/xm/4.5/im ratios at radial distances greater than 
20" (Fig. [6}. We also note that the 5.8/xm/4.5^m and 
8.0/im/4.5/xm ratios take minimum values at a radial dis- 
tance ~15 /; , i.e., the position of the bright [N 11] and H2 
ring. 
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Figure 5. (left) Colour-composite NOT [O in] (blue), Ha (green), and [N n] (red), and (right) Spitzer IRAC 4.5 /zm (blue), 5.8 /im 
(green), and 8.0 /zm (red) pictures of NGC 6369. We note that the stretch of the NOT optical images has been chosen to highlight the 
faint emission from the knots and condensations outside the ionised nebular shell. We also note that the north-south band of 8.0 /zm 
emission is an instrumental artifact. 



3.2 The extended emission 

The HST WF3 images (Fig. [T]) hint the presence of irregular 
condensations outside the envelope of NGC 6369, towards 
the east and close to the edge of the WF3 FoV. The NOT 
[N n] image (Figure 0-Ze/fc), which covers a greater FoV than 
the HST images, and very especially the WHT H2 image 
(Figure. reveal patches of diffuse emission and blobs, 

as marked by arrows in the H2 image, well beyond the main 
ionised nebula, up to ~70" from the central star in [N 11] 
and ^80" in H 2 . These [N 11] and H 2 blobs are distributed 
irregularly, mostly along the east- west direction, with some 
hints of point-symmetry. These blobs are also detected in 
the Spitzer IRAC 8 fim image (Figure fright), especially 
those that are well beyond the 8 /im bright envelope. 



4 MID-IR SPECTROSCOPY 

The ISO spectrum of NGC 6369 in the 2.4-4.1 /zm spectral 
range shown in Figure [8] reveals the prevalence of emission 
lines, with Bra being the most intense and relatively weak 
or negligible emission of the 3.3 fim PAH band. There is 
also evidence for the H 2 1-0 Q(3) A2.424 fim and H 2 1-0 
Q(7) A2.500 /im lines, which is consistent with the detection 
of H 2 emission in our WHT L IRIS images and at 17.036 
/im dPottasch &; Bernard- S alas] l2008h . Although the spec- 
trum remains relatively noisy, and some possible transitions 
are unidentified, we see evidence for Pf£ and Pf7, as well as 
a couple of Humphrey transitions. 

The longer wavelength ISO spectra are relatively noisy, 
precluding reliable analysis in the range A A ~4.1-10 /im, 
although the enhancements in emission at 6.2 and 7.8 fim 
suggests that PAH band features are present. The Spitzer 
IRS spectra (Figure [9]) show more convincing evidence for 
these PAH components. The shorter wavelength SL2 spec- 
trum (FigMleft) shows evidence for the H 2 0-0 S(7) A5.511 
fim line, but the spectrum is, otherwise, dominated by four 
PAH bands at 6.2, 6.9, 7.8, and 8.6 /im, and an underly- 
ing continuum which raises towards longer wavelengths and 
that can be attributable to cool dust emission. 

The same situation applies at longer wavelengths as well 
(Fig. fright), where the flux is dominated by PAH band 



emission yet again, with features at 7.8, 8.6, 11.3, 12.7 and 

14.2 jLim contributing a significant fraction of the total emis- 
sion. We also see evidence for strong transitions due to the 
[Ar ill] A8.991 fim and [S iv] A10.51 fim lines, and the pos- 
sible contribution of the Hua A12.372 fim line. 

One final point is related to the relative intensities of the 
PAH bands themselves. The 7.8 /im feature arises primarily 
through C-C stretching modes, whilst the 11.3 /im band 
is a consequence of C-H out-of-plane bending modes. This 
would normally lead to I(7.8/xm)/I(11.3/im)«1.3 for neutral 
PAHs, but radiation fields in PNe are capable of appreciably 
ionising the PAH molecules, reducing the strength of the 

11.3 /im feature and leading to I(7 1 8m n)/I(11.3/im)^12 for 
fully ionised PAHs (iLi fe Dra ine 2001). It is thus puzzling 
that the I(7.8/xm)/I(11.3/xm) ratio appears significantly less 
than unity (Fig. fright) at any location of the nebula. 

4.1 Interpreting the mid-IR images of NGC 6369 

The mid-IR spectroscopy of NGC 6369 described above al- 
lows us to discuss the nature of the emission detected in 
the IRAC bands. The ISO data show convincingly that 
the emission in the IRAC 3.6 and 4.5 /tm bands is largely 
dominated by H 1 lines, with the Bra A4.052 /mi line in 
the 4.5 fim IRAC band being the brightest line. The con- 
tribution to these bands from Bremsstrahlung and the 3.3 
/tm PAH band is very small or negligible. The contribution 
from other ionic lines in the 4.1-5.1 /tm spectral range (e.g., 
[Mg iv] A4.486 /im, [Ar vi] A 4.529 /zm, and [K in] A4.618 
/xm, iBernard Salas et~al1l200lh is also possible, but the ISO 
spectrum has too low S/N and the Spitzer IRS SL2 spec- 
trum missed this spectral range. Meanwhile, the emission 
in the 5.8 and 8.0 /im IRAC channels is dominated by the 
6.2, 6.9, 7.8, and 8.6 /im PAH bands, the [Ar 111] A8.991 /im 
emission line, and an underlying continuum from cool dust 
which raises towards longer wavelengths. 

The close resemblance between the images and surface 
brightness profiles in the 3.6 and 4.5 /im IRAC channels is 
due to the prominence of ionic and H 1 transitions in these 
two bands. The higher relative emission in the 5.8 and 8 
/im bands at large radial distances implies that molecular 
material and cool dust, responsible of PAH bands and con- 
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Figure 6. MIR IRAC profiles of NGC 6369 along approximate E-W (left) and N-S (right) axes through its central star as indicated in 
the inserted images. The lower panel shows flux ratios for the central region of the source, where 4.5 /im emission has a reasonably high 
S/N. 



tinuum emission, are important at the eastern and western 
extensions of the inner shell and at the envelope, whilst ionic 
lines dominates the emission from the bright inner shell. The 
relative increase of 5.8 and 8.0 /mi, [N n], and H2 emission 
at the interior of the bright inner shell most likely reveal 
an opening of the inner shell at its polar caps, so that the 
emission detected at this region comes from regions far from 
the nebular center, above and below the bright inner shell. 

Finally, we noted in fc|3.2l that the blobs and knots of dif- 
fuse emission outside the main nebular shells are detected 
in the optical [N 11], near-IR H2, and IRAC 8 /mi images. 
Since these blobs present both emission from ionised (N + ) 
and molecular (H2) material, and the IRAC 8 /mi band of 
the main nebula is dominated by PAH bands and dust con- 
tinuum, the origin of the emission in the IRAC 8 /mi band 
may be mult i- fold. The [Ar ill] A8.991 /im line is the bright- 
est ionic line in this band, but the ionization potential of 
Ar+ (IP=27.6 eV) is higher than this of N (IP=14.5 eV) 
and thus Ar ++ may be not present in these regions. Indeed, 
the lack of [O ill] emission, whose IP(0 + ) is 35.1 eV, in- 
dicates that only low- excitation species are present in these 
regions. Since no emission from these blobs is detected in the 
IRAC 5.8 /im image, it seems unlikely that PAH bands or 
dust continuum produce the emission detected in the IRAC 
8 /xm image, although we concede that the levels of dust con- 
tinuum and PAH emission in the IRAC 5.8 /im may be lower 



than in the IRAC 8 /im band. Given the previous difficulties 
for ionic, PAH or continuum emission and that near-IR H2 
1—0 S(l) emission is detected, we favor the H2 0—0 S(4) 
A8.0251 /im emission line as the responsible for the emission 
from these blobs in the IRAC 8 /im image. 



5 DISCUSSION 

5.1 Spatio-kinematical modelling of the inner 
shell of NGC 6369 

Figure displays the individual Ha and [N 11] position- 
velocity (PV) maps of NGC 6369 obtained from the five 
slits positions illustrated in Fig. [H Overall, the Ha and 
[N 11] lines show similar kinematical features, but these are 
best seen in [N 11] given that its thermal broadening is 
smaller than that of the Ha line. The radial velocity of the 
source with respect to the Local Standard of Rest was deter- 
mined to have a value Vlsr = —89.5 =b l.Okms -1 , based on 
the mean expansion velocity at the position of the central 
star. This result is in excellent agree ment with the value o f 
Vlsr — — 90±6.5kms _1 reported bv lSchneider et aTl ([T983). 

The line shapes of the echellogram at slit position #2, 
which cuts the nebula along the eas t and west extensions, i s 
closely similar to that presented bv lSteffen fe Loped (2006): 
the east and west extensions are detected as red- and blue- 
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Figure 7. (top) Ha and (bottom) [N n] A6583 position-velocity (PV) maps of NGC 6369 for slits positions #1 to #5. The slit orientation 
is labeled on the upper panels. The relative position is measured with respect to the horizontal dashed lines that mark the midpoint of 
the most intense central emission for slits #1, #3, #4, and #5, and the central star for slit #2. The vertical dashed lines correspond to 
the systemic radial velocity (Vlsr = —89.5 db lkms -1 ). 



shifted loops, respectively, and their shape is typical of the 
emission in bipolar outflows. The structure of the Ha line 
may give the false impression that the inner shell is a closed 
ellipsoid, but the sharper view offered by the [N n] line re- 
veals that the loops of the east and west extensions emanate 
from gaps at the tips of the inner shell. 

The eastern extension is mapped in more detail by the 
slits at positions #3, #4, and #5. In the echellograms at po- 
sitions #4 and #5, this extension is detected as a red-shifted 
loop, as expected for a blister- like structure. However, the 
echellogram at position #2 reveals a waist and two open 
loops that are indicative of a more complex structure. 

Previous three-d i mensional modeling of NGC 6369 
(|Monteiro et all 12004 ISteffen fc Loped l2006h concur that 
this nebula has a bipolar stru cture, although the d etails 
make them differ sie: nificantly. iMonteiro et all |2004) ob- 
tained spatially-resolved information of the physical condi- 
tions within the nebula that, in conjunction with a photo- 
ionization model, allowed them to synthesize narrow-band 
images. In order to match the morphology observed in sev- 
eral emission lines, they adopted for NGC 6369 a clumpy 
hourglass structure with a thick equatorial belt and two 
symmetrically located caps of emission whose orientation 



is tilted by 30° from the symmetry axis o f the main 
nebula. On the other hand, ISteffen &; Loped (|2006h used 
a single high-dispersion long-slit spectrum of NGC 6369 
obtained along its major axis to illustrate the capabili- 
ties of "SHAPE", an interactive 3-D modeling tool. Us- 
ing a Hubble-type (v(xr) expansion la w for NGC 6369, as 
it is typically assumed in PNe (e.g ., Solf &; Ulr ich 1985; 



[ Guerrero, Manchado, &; Serra-Ricartl Il99 6: Sabb adin et alj 
[2000; ICorradill2004l ) , they concluded that there is no equato- 
rial waist in NGC 6369, but its three-dimensional structure 
can rather be described by a partial spheroidal shell with 
opened tips at its major axis from where two fainter blisters 
protrude as the east end west extensions. 

Our more detailed spatial c overage of the ki n emati cs 
of NGC 6369 allows us to refine ISteffen fe Loped (|2006l Vs 
model and to investigate the spatio-kinematical structure 
of the nebular envelope. Using the latest version of the 
code SHAPE (version 10: ISteffen et alJl201lh . we have con- 
structed a spatio-kinematical model for NGC 6369 which fits 
its narrow-band images and PV maps. The quality of the fit 
can be judged in Figure [TUl that presents a comparison be- 
tween the synthetic morphology and kinematics from the 
model and those observed in NGC 6369. The synthetic im- 
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Figure 8. 7*90 spectrum of NGC 6369 extending between 2.4 
and 4.1 jum, together with the ranges corresponding to the 3.6 /im 
(solid profile) and 4.5 fim (dashed profile) Spitzer IRAC channels. 



age of NGC 6369 (right image of the upper panel in Fig. fTQ|) 
is well matched with its HST Ha image (left image). Sim- 
ilarly, the kinematics of NGC 6369 are fairly reproduced as 
shown in the lower panel of Fig. [10] which presents the PV 
for slits #2 to #5 accompanied by the corresponding PV 
map extracted from the model. Our best fit model, illus- 
trated in Figure [Til is composed of: 

• an equatorial barrel with expansion velocity ~55 km s _1 
and whose axis of symmetry is inclined by ^50° with respect 
to the line of sight (RLoS), 

• a southern lobe approaching to us with an inclination an- 
gle ^40° (RLoS) and polar expansion ~60 km s _1 , and 

• two northern lobes with inclination angles ~115° and 
~127° receding from us with polar expansion velocities ~77 
and ~84 km s _1 , respectively. 

In this way, the physical structure of NGC 6369 is remi- 
niscent of other barrel-shaped PNe with po lar extensions 
that result in a spindle shape (e.g., NGC 3918. ICorradi et ah! 
Il999h . 

Accordin g to this model and fo r the distance of 1,550 pc 
proposed bv lMonteiro et alj (|2004r h the kinematical ages of 
the barrel and bipolar lobes are ~2,000 yrs and ~3,200 yrs, 
respectively. The kinematical age of the inner shell is thus 
smaller than that of the bipolar lobes, most likely indicating 
important dynamical effects that have made the kinemati- 
cal and real ages to differ. These dynamical effects may also 
be revealed by the distorted eastern structure whose com- 
plex shape would arise from the interaction of an expanding, 
regular lobe with extended and irregular material along its 
direction. Alternatively, this bifurcated structure may arise 
from intrinsic multipol ar collimated outflows as has been 
claimed in other PNe (Ma nchado, Stanghellini, &; Guerrero! 
ll996l : ISahaill2000h . 



5.2 Kinematics of the envelope 

The emission from the outer envelope of NGC 6369 is regis- 
tered by the slits at positions #1, #3, #4, and #5, although 
the echellograms at positions #3 and #4 only cover a small 
fraction of this shell. In the echellograms at positions #1 
and #5, the Hoj and [N n] lines are extremely narrow, with 



FWHM ~ 32kms _1 and ~ 20 km s" 1 , respectively. We 
note, however, that the velocity centroid along these lines is 
not constant: the echellogram at position #1 shows a weavy 
pattern (best seen in the [N n] line), whereas at position 
#5 displays an arc-like shape. We measure velocity differ- 
ences up to 17 km s _1 . We also note that emission is patchy, 
especially in the [N n] line. 

Clearly the envelope is not inert, but the nar- 
row and weavy /arced line shape does not comply with 
the expectations for an expanding shel l, even if it 
had the low expansion velocity (jMeaburn et aD [l996; 
iGuerrero. Villaver. fc Manchadd 1998) presumed for the 
slow AGB wind. The knotty and fragmented appearance 
of the envelope of NGC 6369 indicate that the kinematics 
revealed by the echellograms is governed by the emission 
from individual knots, rather than by the emission from a 
continuous shell. This situation is reminiscent of many other 
PNe where bright low-ionization knots are embedded within 
a faint envelope around a bright inner shell, e.g., NGC 239 2 
and NGC 7354 (|0'Dell et alj 120021 : IContreras et alJl201Ch . 
As for these two PNe, the relatively low velocity differences 
measured for the knots of the envelope of NGC 6369 would 
suggest that they lay at different latitudes with respect to 
the symmetry axis of the barrel, but mostly on a weavy disk- 
like or flattened structure, as for NGC 2392 and NGC 7354. 



5.3 Extension of the photo-dissociation region 

The relative sizes of the regions dominated by ionized ma- 
terial and molecules and dust are illustrated in Figure [12] 
that displays the normalized surface brightness profiles in 
the HST Ra and [N n], WHT LIRIS H 2 , and Spitzer IRAC 
4.5 and 8.0 fim. images. The comparison between these pro- 
files confirms that the surface brightness profiles of the inner 
shell of NGC 6369 of the optical Ha emission line and mid- 
IR IRAC 4.5 and 8.0 fim bands are very similar in shape 
and extent. The mid-IR and optical profiles start to depart 
from each other at the edge of the inner shell: the emission 
from Ha and [N n] drops considerably in an annular region 
between 20" and 25" (i.e., the inner regions of the enve- 
lope), while the emission in the IRAC 4.5 and 8.0 fim bands 
stays higher and the H2 profile shows increased emission. 
On the contrary, on the 25 // -35 // outer rim of the envelope, 
the emission in the Ha and [N 11] lines is enhanced, while 
the emission in the IRAC 4.5 and 8.0 fim bands flattens and 
then decays smoothly. There is no H2 emission in this range 
of radial distances. 

The comparison of the relative shapes and extents of 
these different surface brightness profiles is relevant to the 
understanding of the physical structure of NGC 6369. Its 
inner shell is mostly ionized and thus the emission in the 
IRAC 4.5 and 8.0 fim bands is dominated by ionic emission. 
The shell is enclosed by a series of "tadpole" or "cometary" 
knots whose prominent low- excitation [N 11] and molecu- 
lar H2 emission at their heads result in bright [N 11] and 
H2 peaks at ~13". The inner shell is surrounded by a re- 
gion dominated by molecular material, the inner rim of 
the envelope between ~17" and ^25". This region, outside 
the bright inner shell and beyond the annular region de- 
fined by the tails of the "cometary" knots, seems to be a 
real PDR. The relative brightening of the emission in the 
IRAC 4.5 and 8.0 fim bands at these radial distances is very 
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Figure 9. Spitzer IRS spectrum of NGC 6369 in the SL2 short (left) and SL1 long (right) wavelength ranges taken at the four overlapping 
positions shown in the inset: A (dark blue diamonds), B (cyan open dots), C (green triangles), and D (red squares). Positions A and C 
cover the inner shell and north portion of the envelope, whereas positions B and D also cover the inner shell, but the south portion of 
the envelope. Despite the (symmetric) displacement between positions B+D and A+C, the spectra are mostly consistent. The spectral 
range ~ 7.3-7.6 /im corresponds at two overlapping spectral sequences, and results in this regime should be considered as unreliable. 



likely due to molecular H2 emission, although a contribution 
from continuum dust or P AH emissions is also very likely 
as noted in other PNe (e.g., IPhillips fe Ramos-Lariosll2008l . 
2010; Ra mos- Larios &; Phillips 2008). Finally, the outer rim 
of the envelope is again dominated by ionic emission, but 
we note that low- excitation emission dominates. This may 
explain the flattening of the emission in the IRAC 4.5 and 
8.0 /im emission due to an additional contribution of ionic 
emission. 

It is finally worth noting that sources with compara- 
ble He 11 and H 1 Zanstra temperatures are likely to be 
ionisation bound and enveloped by a shell of neutral ma- 
terial. This leads to a correlation between the Zanstra ratio 
Tz(He n)/Tz( H 1), and mid-IR extent outside of the ionised 
shell (see e.g., I Quino-Mendoza, Phillips, &; Ramos-Lariosl 
1201 lh . NGC 6369 seems to be a peculiar case of this correla- 
tion: whereas the emission from molecular material causes a 
larger dimension of its inner shell at mid-IR wavelenghts, the 
presence of ionized material in its singular envelope results 
in very similar final optical and mid-IR dimensions. 



5.4 Shock excitation? 

In 33.11 we noted the large jump in the [N 11] to [O III] ratio 
at the filamentary structures that enclose the eastern and 
western extensions of the inner shell of NGC 6369. Given 
that these structures are at the tips of bubble-like struc- 
tures that may arise from bipolar outflows, it is very likely 
that the unusually high [N n]/[0 in] ratios stem from shock 
interactions of these structures with surrounding material. 
To investigate the structure of these features in more detail, 
we show in Figure [13] a radial profile across a particularly 
bright filament of the eastern extension extracted from the 
HST Ha, [N 11], and [O ill] images and ratio maps. 

The comparison of the surface brightness profiles in- 
dicates that the peak of all these three emission lines are 
almost coincident, but the spatial profile of the [O ill] line is 
noticeable broader. The inner shoulder of this profile, closer 



to the central star (negative offsets in Fig. [T3]) , can be ex- 
plained as the result of a higher photo-excitation. The shoul- 
der of the line farther away from the nebula, however, can- 
not be explained in this way and it very likely reveals the 
presence of a forward shock propagating onto a lower den- 
sity medium. The [O ill] /Ha ratio profile in the lower-panel 
of Fig. [13] allows us to derive a quantitative value for the 
spatial extent of this post-shock region, ~ l'/2. At the dis- 
tance of 1,550 pc, this corresponds to ^3xl0 16 cm, com- 
parable to the post-shoc k cooling widths determined from 
modelling of shocks fe.g. JPopita 1977; S hull fe McKeel 19791 : 
IPittard et ai1l2005h . 



5.5 The external condensations of NGC 6369 

There is a growing number of PNe in the literature that 
show knots or blobs of emission well separated from their 
main nebular shells that do not form a halo. In some cases, 
the kinematics of these kno ts imply high- veloc ity collimated 
outflows, e.g., Flemin g 1 (|Palmer et al.l Il996h or MyCn 18 
(Red man et al.l l2000h , while in some others, e.g., IC4634 
( Guerrero et al.ll2008h , the velocity of the knots is close to 
the systemic velocity, may be due to motions close to the 
plane of the sky. 

As for NGC 6369, we note that one of the blobs located 
outside its inner shell and envelope is detected in the Ha 
and [N 11] emission lines at the slit position #1 (Fig. [4]). 
This blob, at an offset of ~70" from NGC 6369 central star, 
has a relatively small velocity shift, ^+17 km s _1 , with re- 
spect to the systemic velocity, and its velocity width is small, 
FWHM ~ 26kms _1 in Ha and ~ 21 km s" 1 in [N 11]. 

The kinematical properties of this outer condensation 
make NGC 6369 closer to the case of the lower velocity knots 
of IC4634 than to a case of high- velocity collimated out- 
flows. Meanwhile, the velocity shift, although small, would 
preclude that these outer condensations are part of an in- 
complete halo (as, e.g., the case of IC 4593, iCorradi et al.1 
Il997f ). The origin of these external condensations is uncer- 
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Figure 10. (top) HST Ha (left) and synthetic Ha (right) images of NGC6369. (bottom) PV maps from the model compared to the 
observed PV maps. The synthetic PV maps have been displaced towards higher velocities (to the left in the PV maps) to allow a fair 
comparison. 



tain, but they may probe earlier, point-symmetric ejections, 
or arise from shocks associated with the east- west that excite 
material surrounding the nebula. 

5.6 Physical structure of NGC 6369 

The spatio-kinematical model of the inner shell of 
NGC 6369, the kinematics of its envelope, and the intrincate 
distributions of ionized, molecular, and mid-IR emission can 
be used to obtain a coherent view of its structure. The in- 
ner shell and extensions of NGC 6369 can be described as a 
tilted equatorial barrel with bipolar extensions. 

This inner shell is encircled by a ring of [N n] -bright 
"comet ary" knots. The envelope is formed by an ensemble of 
discrete low-ionization and molecular blobs and knots prob- 
ably embedded within a much fainter, fully ionized shell. 
These knots seem to be distributed at slightly different lati- 
tudes with respect to the barrel symmetry axis, but mostly 
on a weavy disk-like or flattened structure. 

Near- and mid-IR H2 emission encompasses the inner 
shell, revealing a real PDR that seems to be inside the enve- 




Figure 11. Main structures of NGC 6369 represented in a 3D 
model. The left panel corresponds to the inclination and PA pro- 
posed. The right panel is the same structure but rotated to show 
a complementary view of the model. 



lope, a shell of ionized gas. It is unclear whether this PDR 
exists within the envelope, or whether this is a projection 
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Figure 12. Surface brightness profiles of the Ha (black), [N il] 
(blue), WHT H 2 (cyan), IRAC 4.5 /xm (green), and IRAC 8.0 /xm 
(red) along PA 330° through the central star of NGC 6369. The 
profiles have been normalized to the emission peak at the bright 
inner shell. The location of the inner shell end envelope edges are 
marked by vertical dotted lines. 



effect, so that the PDR is located at high nebular latitudes, 
above and below the envelope that would be a rather flat 
structure only present at low latitudes. 



6 CONCLUSIONS 

Using narrow-band optical and near-IR images, broad-band 
mid-IR images, optical long-slit echelle spectra and mid-IR 
spectra, we have investigated the physical structure of the 
planetary nebula NGC 6369. The inner shell and its east- 
west extensions can be modeled in terms of a tilted barrel of 
ionized gas and bipolar outflows at its tips, in broad agree- 
ment with previous observations of the source. The nebula 
is optically thick to H ionizing radiation, thus allowing the 
existence of a PDR that is confirmed by its H2 emission and 
suggested by an additional contribution of dust continuum 
and PAH band emission to the IRAC images in the 5.8 and 
8 fim bands. The envelope does not seem to be a real shell 
enclosing the inner shell, but rather a flattened structure 
at the equatorial regions of NGC 6369. Molecular and low- 
excitation emission of [N 11] is also detected in small irregular 
condensations well outside the main nebular shells. Finally, 
we note that HST images show the presence of "cometary" 
knots at the rim of the inner shell whose appearance is sim- 
ilar to those of features in NGC 2392 and NGC 7354, among 
other PNe. 
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Figure 13. (top) [N 11], Ha, and [O ill] surface brightness profiles 
through a likely shock structure in the eastern extension of the 
inner shell of NGC 6369. Negative offsets are closer to the central 
star, while positive offsets corresponds to regions farther away 
from the main nebula, (bottom) [N 11] /Ha and [O ill] /Ha ratio 
profiles. 



SPM, wich is operated by the Instituto de Astrononria of the 
Universidad Nacional Autonoma de Mexico. 

The William Herschel Telescope is operated on the is- 
land of La Palma by the Isaac Newton Group in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto 
de Astroffsica de Canarias. 

Based on observations made with the Nordic Optical 
Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astroffsica de Canarias. 

This research has made use of the NASA/ IPAC In- 
frared Science Archive, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Admin- 
istration. 

We would like to dedicate this paper in memory of our 
colleague and friend, Prof. John Peter Phillips, who recently 
passed away. 



ACKNOWLEDGMENTS 

GRL acknowledges support from CONACyT and PROMEP 
(Mexico). MAG and GRL are partially funded by grant 
AYA2008-01934 of the Spanish Ministerio de Ciencia e In- 
novation (MICINN). RV, MAG and GRL thank support by 
grant IN109509 (PAPIIT-DGAPA-UNAM) . We also would 
like to thank to the referee R. Corradi who made very ex- 
cellent comments for the improvement of the paper. 

Based on observations made with the Observatorio As- 
tronomico Nacional at the Sierra de San Pedro Martir, OAN- 



REFERENCES 

Acosta Pulido J. A., et al, 2003, INGN, 7, 15 
Aller L.H., & Keyes CD., 1987, ApJS, 65, 405 
Bernard Salas J., Pottasch S. R., Beintema D. A., Wesselius 

P. R, 2001, A&A, 367, 949 
Cohen M., & Barlow M. J., 2005, MNRAS, 362, 1199 
Contreras M. E., Vazquez R., Miranda L. F., Olgum L., 

Zavala S., & Ayala S., 2010, AJ, 139, 1426 
Corradi R. L. M., Guerrero M., Manchado A., Mampaso 

A., 1997, New A, 2, 461 



14 G. Ramos-Larios, M.A. Guerrero, R. Vazquez & J. P. Phillips 



Corradi, R. L. M., Perinotto, M., Villaver, E., Mampaso, 

A., k Goncalves, D. R. 1999, ApJ, 523, 721 
Corradi R. L. M., 2004, ASPC, 313, 148 
de Graauw T., et al., 1996, A&A, 315, L49 
Dopita M.A., 1977, ApJS, 33, 437 
Fazio G., et al., 2004, ApJS, 154, 10 

Guerrero M. A., Manchado A., k Serra-Ricart M., 1996, 
ApJ, 456, 651 

Guerrero M. A., Villaver E., k Manchado A., 1998, ApJ, 
507, 889 

Guerrero M. A., et al., 2008, ApJ, 683, 272 
Hajian A. R., Balick B., Terzian Y., k Perinotto M., 1997, 
ApJ, 487, 304 

Hony S., Waters L. B. F. M., k Tielens A. G. G. M., 2001, 

A&A, 378, L41 
Houck J.R., et al., 2004, ApJS, 154, 18 
Hrivnak B. J., Volk K., k Kwok S., 2009, ApJ, 694, 1147 
Kwok S., Volk K., k Hrivnak B. J., 1999, IAUS, 191, 297 
Li A., k Draine B.T., 2001, ApJ, 554, 778 
Manchado A., Stanghellini L., k Guerrero M.A., 1996, 

ApJ, 466, L95 

Meaburn J., Lopez J. A., Bryce M., k Mellema G., 1996, 

A&A, 307, 579 
Meaburn J., Lopez J. A., Gutierrez L., Quiroz F., Murillo 

J. M., Valdez J., Pedrayez M., 2003, RMxAA, 39, 185 
Monteiro H., Schwarz H.E., Gruenwald R., k Heathcote S., 

2004, ApJ, 609, 194 
O'Dell C. R., Balick B., Hajian A. R., Henney W. J., Burk- 

ert A., 2002, AJ, 123, 3329 
O'Dell C. R., Weiner L. D., k Chu Y.-H., 1990, ApJ, 362, 

226 

Palmer J. W., Lopez J. A., Meaburn J., k Lloyd H. M., 

1996, A&A, 307, 225 
Phillips J. P., 2003, MNRAS, 344, 501 
Phillips J.P., k Ramos-Larios G., 2008, MNRAS, 383, 1029 
Phillips J.P., Ramos-Larios G., 2010, MNRAS, 405, 2179 
Pittard J.M., Dobson M.S., Durisen R.H., Dyson J.E., 

Hartquist T.W., k O'Brien J.T., 2005, A&A, 438, 11 
Pottasch S.R., k Bernard-Salas J., 2008, A&A, 490, 715 
Quino-Mendoza J. A., Phillips J. P., Ramos-Larios G., 

2011, RMA&A, 47, 31 
Ramos-Larios G., k Phillips J.P., 2008, MNRAS, 390, 1014 
Redman M. P., O'Connor J. A., Holloway A. J., Bryce M., 

k Meaburn J., 2000, MNRAS, 312, L23 
Sabbadin F., Cappellero E., Benetti S., Turatto M., k 

Zanin C. 2000, A&A, 355, 688 
Sahai R., 2000, ApJ, 537, L43 

Schneider S.E., Terzian Y., Purgathofer A., k Perinotto 

M., 1983, ApJS, 52, 399 
Schwarz H. E., Corradi R. L. M., Melnick J., 1992, A&AS, 

96, 23 

Shull J.M., k McKee C.F., 1979, ApJ 227, 131 

Solf J., k Ulrich H., 1985, A&A, 148, 274 

StefTen, W., Koning, N., Wenger, S., Morisset, C, k Mag- 

nor, M. 2011, IEEE Transactions on Visualization and 

Computer Graphics, 17, 454 
StefTen W., k Lopez J.A., 2006, RMA&A, 42, 99 
Van Winckel H., k Reyniers M., 2000, A&A, 354, 135 
Zhang K., Jiang B. W., k Li A., 2009, MNRAS, 396, 1247 
Zuckerman B., k Aller L.H., 1986, ApJ, 301, 772 



This paper has been typeset from a TgX/ E^T^X file prepared 
by the author. 



